Disordered Potts model on the diamond hierarchical lattice: 
Numerically exact treatment in the large-g limit 
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We consider the critical behavior of the random g-state Potts model in the large-g limit with 
different types of disorder leading to either the nonfrustrated random ferromagnet regime or the 
frustrated spin-glass regime. The model is studied on the diamond hierarchical lattice for which 
the Migdal-Kadanoff real-space renormalization is exact. It is shown to have a ferromagnetic and 
a paramagnetic phase and the phase transition is controlled by four different fixed points. The 
state of the system is characterized by the distribution of the interface free energy P{P) which is 
shown to satisfy different integral equations at the fixed points. By numerical integration we have 
obtained the corresponding stable laws of nonlinear combination of random numbers and obtained 
numerically exact values for the critical exponents. 

PACS numbers: 05.50.-|-q,64.60.F-,75.50.Lk 



I. INTRODUCTION 

Despite continuous efforts, several properties of many 
body systems in the presence of quenched disorder and 
frustration are still not well understood. Notoriously dif- 
ficult systems are spin glasses, in particular, in finite di- 
mensions and with finite-range interactionsi^ One basic 
problem of our understanding in this field of research is 
the lack of exact solutions for nontrivial models. Ex- 
ceptions in this respect are such systems in which dis- 
order fluctuations are fully dominant and the properties 
of the system are governed by an infinite-disorder fixed 
point.^ This happens, among others for random quan- 
tum spin chains and for a class of stochastic models with 
quenched disorder. If, however, the properties of the sys- 
tem are controlled by a conventional (finite-disorder) ran- 
dom fixed point, such as for classical spin glasses, the 
exact results are scarce. 

In these cases, in order to obtain more accurate infor- 
mation, one often considers hierarchical lattices'^ in which 
the Migdal-Kadanoff renormalization^ can be performed 
exactly. For simple models, such as for directed polymers, 
one can notice a simple, although nonlinear relation be- 
tween the original and the transformed random energies 
and one can derive closed integral equations between the 
corresponding distribution functions, which can then be 
studied by various methods. ^'^ However for spin mod- 
els, such as the random-bond Ising model, it is generally 
not possible to write renormalization equations directly 
for the distribution functions. In these cases one can 
treat numerically large finite samples exactly and average 
the obtained results over quenched disorder i^'^ i-'^^'^^i^^ 
Although bringing very accurate numerical results, this 
type of treatment is still not exact and, as usually hap- 
pens in random systems, the source of inaccuracy comes 
from the averaging process over quenched disorder. In 



such calculations one cannot, for example, decide about 
the universality of the fixed points, i.e., if the critical sin- 
gularities are independent or not of the form of the initial 
distribution of the disorder. 

In this paper we consider such a spin model, the g-state 
random Potts model^'^ in the large-q limit, for which the 
Migdal-Kadanoff renormalization leads to closed integral 
equations in terms of the distribution functions. In this 
respect our results are comparable with those obtained 
for directed polymers on hierarchical lattices.— The q- 
state Potts model, as it is well known, is equivalent for 
q = 2 to the Ising model. On a hypercubic lattice for 
sufficiently large q the transition turns to first order, 
whereas on a hierarchical lattice the transition for any 
finite q stays of second order, but the critical exponents 
are q dependent. 

Properties of the Potts model with random ferromag- 
netic and antiferromagnetic couplings (± J model) have 
been studied numerically to some extent. On the cu- 
bic lattice a spin-glass (SG) phase has been identified 
for q = i (Ref. 15) but the SG phase is absent for 
large enough value of q, as observed for q = 10.— On 
the square lattice, the dimension being below the crit- 
ical one, dc > 2, there is no SG phase, and the phase 
diagram consists only of the ferromagnetic and the para- 
magnetic phases. The transition between these phases 
has been studied numerically for q = 3 where it is found 
to be controlled by four different fixed points-^ - (see also 
in Ref. ^*). The model has been also considered on the 
diamond hierarchical lattice with random ferromagnetic 
couphngs.ia<20^^ 

Here we extend these studies by including frustration, 
too, and investigate the large-g limit. We note that 
in the disordered case the large-g limit is generally not 
singular the critical behavior is qualitatively simi- 
lar to that of finite q and also the critical exponents are 
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in terms of the Potts-variables cr^ = 1,2, ... ,q, associ- 
ated with the sites indexed by i of the lattice. Here 
the summation runs over nearest-neighbor pairs and the 
couplings Jij are independent and identically distributed 
random numbers, which can be either positive or neg- 
ative. We consider the lavge-q limit, in which case it is 
convenient to rescale the temperature^^ asT' = T In q, so 
that = q^ where /3' = 1/T' {ks = 1) and consider the 
high-temperature series expansion of the model, in which 
the partition function Z is dominated by one diagram: 



subleading terms 



(2) 



FIG. 1: (Color online) First steps of the construction of the 
hierarchical diamond lattice with 6 = 2. 



smoothly saturating as g ^ cxd^ We mention that it 
has been recently suggested that, in the large-g limit, the 
Potts model is a plausible model for supercooled Hquids. 
This is indeed the case within the mean-field approach^^ 
but it seems to be completely different for the nearest- 
neighbor model] ^^'^^ 

To consider the large-g limit of the model leads to tech- 
nical simplifications, at least for random ferromagnetic 
couplings. The high-temperature series expansion of the 
model is dominated by a single diagram, the properties 
of which have been studied in 2c? and 3c? by combinato- 
rial optimization methodsi^ This type of simplification 
is valid for the hierarchical lattice, too, but for this lat- 
tice the exact renormalization holds for antiferromagnetic 
couplings, too. As we will see the renormalized param- 
eters in this limit are expressed as simple but nonlinear 
combinations of the original parameters, such as the in- 
terface free energy. This makes it possible to write in- 
tegral equations for the distribution functions which are 
studied by various methods. 

The structure of the paper is the following: the q- 
state Potts model and its Migdal-Kadanoff renormaliza- 
tion in the large-g limit is presented in Sec. |TT1 The 
phase diagram of the model is calculated by the numer- 
ical pool method in Sec. IIIIl whereas the properties of 
the fixed points are obtained through numerical integra- 
tion in Sec. II VI Our results are discussed in Sec. IVl and 
details of the derivation of the integral equations at the 
fixed points are given in the Appendices. 



II. MODEL AND RENORMALIZATION 



which is related to the free energy through (f> = —j3'F. 

In this paper the Potts model is considered on the di- 
amond hierarchical lattice, which is constructed recur- 
sively from a single link corresponding to the generation 
n = 0. (see Fig. [T]). The generation n = 1 consists of b 
branches in parallel, each branch containing two bonds in 
series. The next generation n = 2 is obtained by applying 
the same transformation to each bond of the generation 
n = 1. At generation n, the length L„ measured by the 
number of bonds between the two extreme sites A and B 
is Ln — 2", and the total number of bonds is 



with c?e//(6) = (3) 



where deff(b) represents some effective dimensionality. 

In the following we consider two different boundary 
conditions and denote the partition function as Z^'^ and 
-E,j'^, when the two extreme sites A and B are fixed in 
the same state or in different states, respectively. Their 
ratio is given by: 
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(4) 



where F^^^^" = F^'^ — F^'^ is the interface free energy. 
The ratio Xn obeys the recursion equatio n^^i^^ 



n 



+ (9-1) 
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(5) 



Now we consider the large- g limit and keeping in mind 
that the partition function is given by one dominant 
term, see Eq. ([2]), we obtain the following recursion for 
the scaled interface free energy, /„ — P'F™^"^": 



We consider the g-state Potts model defined by the 
Hamiltonian: 



W = - ^ Jlj5{(7i,(7j) 



In+l =Yj^ 
i=l 



(1) Here the auxiliary function is: 
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with Imax = 

The initial condition is given by 



and = mm(/(i) , 



r(») 



(8) 



where Ji is the value of the ith coupling. 

Note that in the large-g limit the recursion relation in 
Eq. ([5|) are simplified and there is now a direct relation 
in Eqs. ([6|) and ([7]) between the scaled interface free en- 
ergies. This recursion relation involves a somewhat com- 
plicated nonlinear combination of random variables and 
we are interested in the stable law for their distribution 
function. 



III. PHASE DIAGRAM 



For p > all couplings are random ferromagnetic and 
in the limit p — > 1 we have the pure system. For p < 
there are also negative bonds, their fraction is increasing 
with decreasing p. For p = —1 the distribution is sym- 
metric with zero mean and standard deviation 1/VT2- 
In the numerical calculations we have used the so-called 
pool method. Starting with N random variables taken 
from the original distribution in Eq. (|10p . we generate a 
new set of TV variables through renormalization at a fixed 
temperature, T', using Eqs. (O and 0. These are the 
elements of the pool at the first generation, which are 
then used as input for the next renormalization step. We 
check the properties of the pool at each renormalization 
by calculating the distribution of the scaled interface free 
energy, its average and its variance. In practice we have 
used a pool of iV = 5 x 10^ elements and we went up to 
77, ^ 70 — 80 iterations. 



A. Non-random system 



In the pure case the recursion relation simplifies: 





/„+i = { 6(2/„ - 1) 

bin 



if < /„ < 1/2, 
if 1/2 < /„ < 1, 
if /„ > 1. 



(9) 



It has a fixed-point at Ic = 6/(26— 1). The phase transi- 
tion at this point is of first order. The dominant diagram 
for I < Ic consists of isolated points, whereas for I > Ic 
the dominant diagram is fully connected and there is a 
phase coexistence at I = Ic- 



B. Random system: Numerical study 

For the random case we consider a diamond lattice 
with branching number 6 = 2, i.e., with an effective di- 
mension dcff — 2. First we perform a numerical inves- 
tigation using a continuous box-like distribution of the 
couplings 



V{J) 
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if < J < T^, 

1— P 1— p' 



otherwise. 



(10) 



1. Phases 

At a given point of the phase diagram, (j),T'), the 
renormalized parameters display two different behav- 
iors, which are governed by two trivial fixed points. In 
the paramagnetic phase the scaled interface free energy 
renormalizes to zero, whereas in the ferromagnetic phase 
it goes to infinity. We note that in the large-g limit there 
is no zero-temperature spin-glass phase, contrary to the 
known results for q = 2 and q = Sfiii^i We analyze the 
properties of these trivial fixed points in details in the 
following section. 

The two phases are separated by a phase transition line 
T^{p) which can be calculated accurately for a given pool 
and its true value can be obtained by averaging over dif- 
ferent pools and taking the N oo limit. The phase dia- 
gram is shown in Fig. [5]in the plane T'/J versus p. In the 
phase diagram one may notice a reentrance in the regime 
p < 0. At low temperatures the system is disordered due 
to frustration. In the intermediate temperature regime 
there is an order-through-disorder phenomenon and at 
high temperature the system is again disordered due to 
thermal fluctuations. Similar features has been analysed 
before for the q — 3 modeli^ 



with p < 1 and such that the mean value of the couplings 
is given by: 



2. Phase-transition lines 
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(11) 



The properties of the phase transition are different 
when the disordering effect is dominantly of thermal ori- 
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gin (which happens in the high-temperature part of the 
transition Hne) or when it is dominantly due to frustra- 
tion (which happens in the low-temperature part of the 
transition hne). At the boundary of the ferromagnetic 
part of the phase diagram, the pure-system fixed point 
at p = 1 is unstable in the presence of any amount of (fer- 
romagnetic) disorder and the transition is controlled by 
a new fixed point, the random- ferromagnet (RF) fixed 
point. According to the phase diagram in Fig. [2l this 
fixed point controls the phase transition even in a part 
of the region where p < 0, up to a point (pmCj T^{pmc))- 
Our numerical studies indicate that at the RF fixed point 
the scaled interface free energy is typically /„ — 0(1). 
Below this temperature, T' < r^(pMc), the phase transi- 
tion is controlled by a zero-temperature (Z) fixed point, 
the properties of which are very similar to that in the 
(7 = 3 model^^, in which it describes the transition be- 
tween the zero temperature spin-glass phase and the 
paramagnetic phase. Our numerical studies show that 
along the green line of Fig. [2] between MC and Z, as well 
as at the Z fixed point the scaled interface free energy 
grows with the size, L, as 

I{L) = L'^^ui , (12) 

where the ui are 0(1) random numbers and the droplet 
exponent is 9z > 0. Finally at (pmCi T|!(pmc)) there is a 
multicritical (MC) fixed point, analogous to the Nishi- 
mori MC point in gauge-invariant systems On the 
phase diagram shown in Fig. [2l the coordinates of the 
fixed points are (-0.212334; 0.353194) for the MC fixed 
point and (—0.173626; 0.0) for the Z fixed point, but we 
did not study the actual position of the RF fixed point. 
In the following section we study the properties of the 
fixed points (except the MC point) through numerical 
solution of integral equations. 



IV. PROPERTIES OF THE FIXED POINTS 

As shown in the numerical study of the previous sec- 
tion, the system has two phases which are controlled by 
two trivial fixed points and the phase transition line is 
controlled by four different nontrivial fixed points. In 
each fixed point there is a characteristic distribution of 
the scaled interface free-energy parameters, P{I), which 
transforms under recursion in Eqs. ([6]) and (O into P'{I). 
The transformation law of the distribution function can 
be written in an explicit form, due to the fact that the 
transformed variables, In+i, are expressed as a sum of 
a few /„, which are independent random numbers dis- 
tributed according to P{I). Thus we are looking for the 
stable law of a nonlinear combination of random num- 
bers, which is different in the different fixed points. 
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FIG. 2: (Color online) Phase diagram of the disordered Potts 
model in the large-g limit on the diamond hierarchical lat- 
tice with b — 2 using the continuous box-like distribution 
of Eq. (|10|) . The critical behavior along the transition line 
between the paramagnetic and the ferromagnetic phases is 
controlled by four different fixed points: i) At p = 1 the pure 
systems fixed point (P), ii) from p = 1 to p = pmc (red line), 
the random- ferromagnet fixed point (RF) , iii) at p = Pmc the 
multicritical fixed point (MC), iv) for pmc < P < Pz (green 
line) the zero temperature fixed point (Z). Note that there is 
a reentrance in the phase diagram. 



A. Paramagnetic phase 

The paramagnetic phase has a trivial limit distribution 
Ppara = ^(I) sincc all parameters renormalize to zero. 



B. Ferromagnetic phase 

In the ferromagnetic phase grows without limit, thus 
asymptotically for each bond Imax > Imin > 1, and in 
Eqs. ^ the last equation holds. Consequently for 6 = 2 
we have: 

/„+i = /,« (mm) + (mm) , (13) 

which leads the following relation for the distribution 
function: 

P'{I) = I dxP2{x)P2{I-x) , (14) 

Jo 

in terms of 

/oo 
dxP{x), I>1. (15) 

Writing I as I + Ii where / is the average value and 
Ii the fluctuating part, Eq. (fT4|) leads to a renormaliza- 
tion relation for the average values, (/)' = 21. Defining 
the probability distributions for the fluctuating part as 
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probability po = 0.1280795 and the probability distribu- 
tion is continuous, but nonanalytic, at / = 1 and 1 = 2. 

In order to calculate the thermal eigenvalue of 
the fixed-point transformation we form the Jacobian 
J{x, y) through functional derivation at the fixed point, 
J{x,y) = SP'{x)/5P{y), and solve the eigenvalue prob- 
lem 
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FIG. 3: (Color online) Probability distribution of the scaled 
interface free energy at the RF fixed point. At / = there 
is a delta peak with strength, po = 0.1280795, the function is 
nonanalytic at J = 1 and 1 = 2. 



P{h) = P{I + h) and P\h) 
the following relation: 



P'(2/ + /i), they satisfy 



P'(/i) = 4 / dxiP{xi)P{h-xi) 



dyiPivi) / dziP(zi) . (16) 

We have studied this equation numerically and found 
that the fixed-point solution satisfies the relation 



P'ih) = \p{h/\) 



(17) 



with A = 1.230091(1). Consequently in the ferromagnetic 
fixed point the average value and the standard deviation 
of the scaled interface free energy are related to the size 
L of the system by: 

- ^d. 



(18) 



Here the dimension of the interface is dg = d^fj — 1 = 
1 and the droplet exponent is 9 = log(A)/ log(2) 
0.298765(1). Note that the value of the droplet expo- 
nent is in agreement with previous numerical studies for 
finite value of It is also identical to the value obtained 
for directed polymers on the same hierarchical latticei^i^ 



C. The RF fixed point 

At the RF fixed point the iterated interface free energy 
scales to the region /„ > 0, therefore in the recursion re- 
lations of Eq. ([7]) the second is irrelevant. The probability 
distribution P{I) satisfies different recursion relations at 
/ = and in the regions 0</<l, l</<2, and / > 2. 
These relations are obtained in Appendix [X] At the fixed 
point the distribution stays invariant, thus P'{I) = P{I)- 
The fixed point distribution obtained by numerical inte- 
gration is shown in Fig. [3] Note that / vanishes with 



dyJ{x,y)My)^K'M^) 



(19) 



In this way we have obtained the leading eigenvalue 
= 1.6994583(1), and checked that < 1 for 

i > 2, e.g., the correction-to-scaling eigenvalue is Aj^ = 
0.6378796(1). 

The thermal eigenvalue of the transformation is given 
by: 



Vt 



RF 



RF 



log A? 
log 2 



0.7650750(1) 



(20) 



from which we deduce the correlation length critical ex- 
ponent J/RF = l/yj^^ = 1.307061(1). This exponent ap- 
pears in the scaling form of the average interface free 
energy I{T',L) in the ordered phase for T' < T'^ along 
the transition line in Fig. O We have 



I{T\L) = 



(21) 



where dg = o^e// — 1 = 1 and the correlation length di- 
verges as iav{T') - (T' - t;^)-"^". Note that the fluc- 
tuation of the interface free energy AI{T' , L) grows with 
the droplet exponent (see Eq. (fT8|)) as: 



AI{T',L) 



L 



(22) 



The associated correlation length ^variT') is proportional 
to ^av{T'), thus there is only one length scale in the prob- 
lem. The situation seems to be different for finite values 
of q, in which case two different correlation length ex- 
ponents are obtained for the average and the fluctuating 
part of the interface free energy, respectively.— 



D. The zero temperature fixed point 

The Z fixed point is at zero temperature and here the 
scaled interface free energy /„ grows to infinity. There- 
fore we use the reduced variable i„ = In/ In in terms of 
which the renormalization group equations in Eq. ([6]) are 
modified as: 



with an+i = In+i/In and 
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1=1 



(23) 
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> 0. 



(24) 



In agreement with the phase-diagram in Fig. [2] these 
equations have two trivial fixed points corresponding to 
the ferromagnetic phase when p > pz with q;„ —^ b and 
to the paramagnetic phase for p < pz, in which case 
ckn — > 0. The nontrivial fixed point at p — pz where 
an az governs the zero temperature transition. 

The probabihty distribution satisfies different re- 
lations for i < 0, i > 0, as well as at i = 0. The cor- 
responding integral equations are presented in Appendix 
IbI We have integrated these equations numerically and 
found that, at the Z fixed point, the probability distribu- 
tion n(i) transforms as 



n'{i) = azU{i) , 



(25) 




with az — 1.10661(1). Consequently, the droplet expo- 
nent at the zero temperature transition (see Eq.([T2])) is 
given by: 



log az 
log 2 



= 0.14615(1) 



(26) 



The probability distribution at the fixed point is shown 
in Fig. m 

We have calculated the Jacobian at this fixed point, 
too. Like in Eq. (|19p . from the corresponding eigenvalue 
problem we have deduced the leading eigenvalue = 
1.49314(1) which gives the thermal exponent 



yf 



.57835(1) 



(27) 



and the correlation-length exponent 
1.72906(1). 

In the ferromagnetic phase the average interface free 
energy scales along the green transition line in Fig. [5] as 



I{T',L) 



[ea.(T')] 



ds-6 



(28) 



with ds = deff — 1 = 1. Similarly, the fluctuation of the 
interface free energy has the scaling form: 



A/(T',L) = 



[e.ar(T')]^-«- 



+ 



(29) 



Note that these scaling relations differ from those at the 
RF fixed point; see Eqs. (HH) and (j^ . respectively. The 
finite-size behavior of these modified forms with ^{L) ^ L 
is in agreement with Eq. (fT^ . The correlation lengths 
^av and ^ 

i,ar diverge as T' — > with the same critical 
exponent v^. 



FIG. 4: (Color online) Probability distribution of the scaled 
interface free energy, i — I /I, at the Z fixed point, as well 
as along the green line of Fig. [2] between MC and Z. At i = 
there is a delta peak with strength, ttq = 0.000158, the 
function is nonanalytic at i = 0. 



E. The MC fixed point 

To complete our study we list here the properties of the 
MC fixed point, which has been obtained numerically us- 
ing the pool method. At the MC point, the interface free 
energy is found to be an 0(1) random number, the dis- 
tribution of which is shown in Fig. [S] Both negative and 
positive couplings are involved in the distribution and 
the average value of the interface free energy as well as 
its fluctuations are larger than at the RF fixed point. In 
the vicinity of the MC point, the scaling of the interface 
free energy takes the form given in Eqs. ([2T|) and (|22p . 
The divergence of the correlation length can be analyzed 
assuming the presence of a multiplicative logarithmic cor- 
rection. The numerical data seem to be consistent with 
the scaling combination 



(30) 



with K « 1.5 and i^mc = 3.61 for the average as well as 
for the standard deviation. 



V. DISCUSSION 

We have studied the random Potts model in the large- 
q limit with such type of disorder which includes both 
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FIG. 5: (Color online) Distribution of the interface free energy 
at the MC point as calculated by the numerical pool method. 
The green (red) symbols are for n — 12 [n — 10) iterations. A 
small fraction no ~ 0.0057 of the samples have zero interface 
free energy. 



the random (nonfrustrated) ferromagnet as well as the 
frustrated spin-glass regime. The model is considered on 
the diamond hierarchical lattice, on which the Migdal- 
Kadanoff renormalization is exact. First we used the nu- 
merical pool method to determine the phase diagram. 
It consists of a paramagnetic phase and a ferromagnetic 
phase, separated by a transition line which is controlled 
by four fixed points. This structure of the phase diagram 
is very similar to that found numerically for the ± J three- 
state Potts model in two dimensions, although in our 
model the zero-temperature spin-glass phase is absent. 

The state of this random system is shown to be 
uniquely determined by the distribution of the scaled 
interface free-energy, P{I). The renormalization group 
transformation is written in the form of integral equa- 
tions for P{I) and we have studied the properties of its 
fixed points by numerical integration. Mathematically 
the above problem is equivalent to find the stable law of 
nonlinear combination of random numbers. 

In the random ferromagnetic phase the probability dis- 
tribution is analogous to that of directed polymerSf^'^ For 
the lattice with h = 2 the droplet exponent is exactly the 
same in the two cases. The nontrivial fixed points govern- 
ing the properties of the phase transition have different 
scaling properties, which arc summarized in Table |T1 

We note that at the RF fixed point the correlation- 
length exponent z^rf corresponds to a mostly thermal 
scaling field, whereas at the ZZfixed point z/z is mainly 
due to a disorder-like scaling field. In the large-g limit, 
the same correlation length exponent governs the scaling 
form of both the average and the fluctuation part of the 
interface free energy. 

Our results are obtained strictly in the large-q limit 
and our analysis is mainly done for a diamond hierar- 
chical lattice with 6 = 2. Here we comment on possible 



TABLE I: Scaling behavior of the average interface free en- 
ergy, 7, and its fluctuations A/, as a function of the linear size, 
L, in the ferromagnetic and paramagnetic phases and at the 
different fixed points: P (pure system); RF (red line between 
P and MC in Figl2|; MC; Z (green hue between MC and Z in 
Fig[2}. The values of the correlation-length critical exponent 
are also indicated. At the P fixed point the transition is of 
first order. 





I 


AJ 




Ferro 




L' 




Para 










P 

RF 
MC 
Z 


6/(26- 1) 
0(1) 
0(1) 




0(1) 
0(1) 


l/d 
1.31 
3.61 
1.73 



extensions of our study in three directions, i) For finite, 
but not too large values of q the phase diagram is sim- 
ilar to ours, for example the SG phase appears only for 
q- values which are not too far from two. The critical 
exponents are expected to have corrections in powers of 
l/q. ii) One can repeat the calculation for larger values 
of 6, which corresponds to a larger effective dimension- 
ality. However the integral equations for the probability 
distribution P{I) become more and more complicated: 
P{I) is a piece-wise function with more and more dif- 
ferent intervals of definition, iii) Finally, the calculation 
can be extended to other quantities such as the inter- 
face energy and the magnetization. For these quantities, 
however, one should work with conditional probabilities, 
which will result in even more complicated equations. 
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APPENDIX A: INTEGRAL EQUATIONS AT THE 
RF FIXED POINT 

The probability distribution at this fixed point is a 
piece-wise function which satisfies different relations in 
the regions 0</<l,l</<2 and I > 2, respectively. 
It has a delta peak at / = with strength po- 

The strength of the delta peak renormalizes as 

Po = Pi , (Al) 

in terms of 

Po = 2po - Po + / dl Pil) [ dx P{x) . (A2) 
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In the region < J < 1 we have the relation 

P'{I) = 2PoPi{I) + [ dx Pi{x)Pi{I ~ x) , (A3) 
Jo 



with 



Pi{I) = 2 / dxP 

J(I+l)/2 



{x)P{I + l-x) 



/oc 
dxP{x), 0<I<1. (A4) 



In the region 1 < / < 2 we have the relation: 



/-I 



(A5) 



P'il) = 2PoP2(/)+ / dx Piix)P2{I ~ X) 

1 

+ / dxPi{x)Pi{I-x) , 

where the auxiliary function P2{I) is defined in Eq. ([15]). 
Finally for / > 2 the renormalization reads as: 

P'{I) = 2PoP2{I)+ I dx Pi{x)P2{I ~ x) 
II 

+ ^ dx P2{X)P2{I ~ x) , (A6) 



APPENDIX B: INTEGRAL EQUATIONS AT THE 
Z FIXED POINT 

The probability distribution at the Z fixed point is a 
piece-wise function which satisfies different relations in 



the regions i < and i > 0, respectively. It has a delta 
peak at « = 0, with strength ttq, which renormalizes as 



I'd = , 



(Bf) 



in terms of 



din{i) 



IIo = 2po - Po ^ 
In the region i < we have the relation 

—U'(i) = 2noni(i)+ / dxIli(x)U2(i - x) 



(B2) 



+ / cia;ni(a;)ni(i - a;) 



(B3) 



in terms of 



/i poo 
dxU{x) + 2n{i) dxU{x), i<0 
-oo <J —i 

(B4) 



and 



U2{i)^2U{i) dxU{x), i>0. (B5) 

J i 

Finally, in the region z > 0, we have the relation 



a 



/oc 
dxn2ix)Ui{i~x) 

+ f dxn2ix)U2ii-x) . (B6) 
Jo 



* Electronic address: igloi@sz£ki.hu| 

^ Electronic address: turban@lpm.uhp-nancy.fr' 
^ For reviews on spin glasses, see K. Binder and A. P. Young, 
Rev. Mod. Phys. 58, 801 (1986); M. Mezard, G. Parisi, and 
M. A. Virasoro, Spin Glass Theory and Beyond (World Sci- 
entific, Singapore, 1987); K. H. Fisher and J. A. Hertz, 
Spin Glasses (Cambridge University Press, Cambridge, 
1991); Spin Glasses and Random Fields, edited by A. P. 
Young (World Scientific, Singapore, 1998). 
^ F. Igloi and C. Monthus, Physics Reports 412, 277-431, 
(2005) 

^ A. N. Berker and S. Ostlund, J. Phys. C 12, 4961 (1979); 
M. Kaufman and R. B. Grifiiths, Phys. Rev. B 24, 496 

(1981) ; R. B. Grifiiths and M. Kaufman, ibid. 26, 5022 

(1982) . 

^ A. A. Migdal, Sov. Phys. JETP 42, 743 1976; L. P. 

Kadanoff, Ann. Phys. N.Y. 100, 359 (1976). 
^ B. Derrida and R. B. Grifiiths, Europhys. Lett. 8, 111 

(1989); J. Cook and B. Derrida, J. Stat. Phys. 57, 89 

(1989). 

^ C. Monthus and T. Garel, J. Stat. Mech.: Theory Exp. 
P01008 (2008) 

A. P. Young and R. B. Stinchcombe, J. Phys. C 9, 4419 
(1976); B. W. Southern and A. P. Young, ibid. 10, 2179 



(1977). 

^ C. Jayaprakash, E. K. Riedel, and M. Wortis, Phys. Rev. 

B 18, 2244 (1978). 
^ S. R. McKay, A. N. Berker, and S. Kirkpatrick, Phys. Rev. 

Lett. 48, 767 (1982); E. J. Hartford and S. R. McKay J. 

Appl. Phys. 70, 6068 (1991). 
^° E. Gardner, J. Phys. France 45, 1755 (1984). 
" A. J. Bray and M. A. Moore, J. Phys. C 17, L463 (1984); J. 

R. Banavar and A. J. Bray, Phys. Rev. B 35, 8888 (1987); 

M. A. Moore, H. Bokil, and B. Drossel, Phys. Rev. Lett. 

81, 4252 (1998); S. Boettcher, Eur. Phys. J. B 33, 439 

(2003). 

M. Nifie and H. J. Hilhorst, Phys. Rev. Lett. 68, 2992 

(1992); M. Ney-Nifie and H. J. Hilhorst, Physica A 193, 

48 (1993; ibid 194, 462 (1993); M. Ney-Nifie, Phys. Rev. 

B 57, 492 (1998); M. J. ThiU and H. J. Hilhorst, J. Phys. 

I 6, 67 (1996). 
13 F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 
1* B. Derrida, C. Itzykson, and J. M. Luck, Commun. Math. 

Phys. 94, 115 (1984). 
1^ L. W. Lee, H. G. Katzgraber, and A. P. Young, Phys. Rev. 

B 74, 104416 (2006). 

C. Brangian, W. Kob, and K. Binder, J. Phys. A 36, 10847 
(2003). 



9 



E. S. Soronscn, M. J. P. Gingras, and D. A. Huse, Euro- 
phys. Lett. 44, 504 (1998). 

J.L. Jacobsen, and M. Picco, Phys. Rev. E 65, 026113 
(2002). 

W. Kinzel and E. Domany, Phys. Rev. B 23, 3421 (1981). 

B. Derrida and E. Gardner, J. Phys. A 17, 3223 (1984); B. 
Derrida, in Critical Phenomena, Random Systems, Gauge 
Theories, edited by K. Ostorwalder and R. Stora, Pro- 
ceedings of the Les Houches Summer School of Theoretical 
Physics, 1984 (North- Holland, Amsterdam, 1986), p. 989. 
D. Andelman and A. N. Berker, Phys. Rev. B 29, 2630 
(1984). 

C. Monthus and T. Garel, Phys. Rev. B 77, 134416 (2008). 
J. L. Cardy, Physica A 263, 215 (1999). 

M. Picco, Phys. Rev. Lett. 79, 2998 (1997); C. Chatelain 
and B. Berche, ibid. 80, 1670 (1998); Phys. Rev. E58, 
R6899 (1998); 60, 3853 (1999); T. Olson and A. P. Young, 
Phys. Rev. B 60, 3428 (1999). 

J.-Ch. Angles d'Auriac and F. Igloi, Phys. Rev. Lett. 90, 



190601 (2003); M.-T. Mcrcaldo, J.-Ch. Angles d'Auriac, 
and F. Igloi, Phys. Rev. E 73, 026126 (2006). 

G. Parisi, in Complex Behavior of Glassy Systems ed M 
Rubi and C Perez- Vicente (Berlin, Springer, 1997) p 111; 
S. Franz and G. Parisi Physica A 261, 317 (1998); M. 
Mezard and G. Parisi, Phys. Rev. Lett. 82, 747 (1999); 
M. Mezard Physica A 265 352 (1999); G. Parisi, Physica 
A 280 115 (2000); X. Xia and P. G. Wolynes Phys. Rev. 
Lett. 86, 5526 (2001). 

I. Morgenstern and K. Binder, Phys. Rev. B22, 288 (1980), 
I. Morgenstern and H. Horner, Phys. Rev. B 25, 504 
(1982), W. L. McMillan, Phys. Rev. B28, 5216 (1983). 
M. J. P. Gingras and E. S. Sorensen, Phys. Rev. B57, 

10264 (1998). 

H. Nishimori, Prog. Theor. Phys. 66, 1169 (1981); M. 
Ohzeki, H. Nishimori, and A. N. Berker, Phys. Rev. E 
77, 061116 (2008). 



